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The exact FZZT brane partition function for topological gravity with matter is computed 
using the dual two-matrix model. We show how the effective theory of open strings on a 
stack of FZZT branes is described by the generalized Kontsevich matrix integral, extending 
the earlier result for pure topological gravity. Using the well-known relation between the 
Kontsevich integral and a certain shift in the closed-string background, we conclude that 
these models exhibit open/closed string duality explicitly. Just as in pure topological 
gravity, the unphysical sheets of the classical FZZT moduli space are eliminated in the 
exact answer. Instead, they contribute small, nonperturbative corrections to the exact 
answer through Stokes' phenomenon. 
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1. Introduction 

Open/closed string duality is one of the most intriguing facets of string theory. It 
states that dynamical processes involving open strings can be formulated strictly in terms 
of closed strings, and vice versa. Often, the open string side of the correspondence can 
be reduced to an ordinary gauge theory. In such cases, the correspondence can provide a 
promising way to address outstanding issues of quantum gravity by reformulating them as 
questions in gauge theory. 

Recently, there has been renewed interest in "minimal string theory," i.e. non-critical 
strings with c < 1. (For a list of references to recent work, see e.g. |1|].) These models are 
dual to certain zero-dimensional gauge theories (i.e. matrix models), providing a tantalizing 
hint that some form of open/closed string duality is at work here. Because these systems 
are integrable, it is natural to expect that c < 1 string theories exhibit open/closed string 
dualities which are explicitly demonstrable. This strongly motivates us to study these 
models further. 

One of the main ingredients in formulating an open/closed string duality are D-branes. 
The D-branes relevant to non-critical strings with c < 1 are called FZZT branes, and they 
were studied in the context of Liouville theory in 0,|3[]. In the recent papers of ||-||, FZZT 
branes were used to provide a geometric interpretation of minimal string theory, as well 
as to probe the structure of its target space. Meanwhile, they were used in |7j to explore 
the relation of minimal strings to topological strings on non-compact Calabi-Yaus with 
B-branes. 

Since our goal here is to understand open/closed string duality in minimal string 
theory, it is natural to consider the dynamics governing the open strings ending on FZZT 
branes. Exactly this question was addressed for the case of the c = —2 model - known 
as pure topological gravity from the work of || - first from the point of view of string 
field theory in || , and later using the double-scaled matrix model in || . Both approaches 
found that the open strings ending on FZZT branes in this background are described by 
the matrix integral of Kontsevich |T(J . By identifying the precise deformation of the closed 
string background associated with the presence of the brane, the authors of |||| showed 
concretely how Kontsevich's formulation of two dimensional gravity can be thought of as 
a kind of open/closed string duality. 

Pure topological gravity is often referred to as the (2, 1) model, by analogy with (p, q) 
minimal string theory. It has natural generalizations, commonly referred to as the (p, 1) 
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models, which describe topological matter coupled to gravity. The appropriate generaliza- 
tions of the Kontsevich matrix integral also exist |ID|jn|] . In this article, we will use the 
double-scaled two-matrix model |H| to relate the two descriptions via open/closed string 
duality. That is, we will show explicitly that the generalized Kontsevich integral describes 
the effective theory of open strings between FZZT branes in the (p, 1) closed-string back- 
ground. Just as for the (2, 1) model, we will see that the generalized Kontsevich integral 
includes nonperturbative effects that drastically modify the topology of the FZZT moduli 
space. 

This paper is organized as follows. We begin in section 2 by formulating the (p, 1) 
minimal string theories as a double-scaling limit of the two matrix model. We compute the 
partition function of a stack of FZZT branes in this background and obtain the generalized 
(matrix) Airy function. We also explore the Stokes and anti-Stokes lines for the FZZT 
partition function. In section 3, we use the correspondence between macroscopic loop 
operators and the local closed string operators to formulate the open/closed string duality 
for the (p, 1) model. Finally, we conclude in section 4 with a discussion of open problems 
and relation to other work. The appendices contain various generalizations and technical 
details. 



2. Double-scaling the (p, 1) models 

2.1. FZZT branes in the two-matrix model 

In this section, we describe the analysis of the FZZT partition function for general 
(p, 1) models using the two-matrix model, extending the results of @. Let us begin by 
considering a general form of the two- matrix model 

Z{g) = J dAdBe-^ v ^ +w(B) - AB) , (2.1) 

where A and B are N x N matrices, g is the coupling constant of the bulk theory, and the 
choice of integration contour depends on the form of the potentials V and W. This model 
can describe (p, q) minimal string theory in the large N limit |12[ , provided we tune the 
potentials V(A) and W(B) to 

^) = si/^fe"-w (22) 
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where we integrate zona small contour around z = 0, and X*(z) and Y*(z) have zeros of 
order q and p, respectively, at z = 1. We find it convenient to take 



XJz) = ^ 



YJz) 



{1-zY 



(2-3) 



z 

Then, for q = 1, the potentials V(a) and W(6) take the simple form 

P k 

V P (°>) = -J2t +Hp > W p( b ) = b - ( 2 - 4 ) 
k=i 

Here we have fixed the integration constant to be 



H P = Y,\ (2.5) 



k=i 



so that V p (l) = 0. This will be convenient for later calculations. One can further simplify 
the potential by shifting A — > A + 1, which leads to 



Z<p,i)(g) = J dAdBe- 1 ^^-^ 



(2.6) 



In the large iV double-scaling limit, this becomes the bulk partition function of the (p, 1) 
model. Two comments on this result are in order: 

1. For p = 2, one can integrate out A to obtain a Gaussian model for B. This case 
was studied in detail in ||. Of course, A cannot be integrated out so easily in gen- 
eral. Nevertheless, the integral (2J5) is essentially trivial for any value of p, since B 



always acts like a Lagrange multiplier constraining A. This will facilitate much of the 
computation reported in this article. 



2. Some care is necessary in order to ensure that the matrix integral in ( |2.6| ) is well 
defined. For even p, this can be accomplished by integrating A and B with respect 
to the measure where r\A and ir/~ 1 B are Hermitian matrices, and if = —1. For odd 
p, one can use the same measure, but one should integrate first over B so that the A 
integral is constrained. 

We will consider exact D-brane observables constructed from the macroscopic loop 

W(y)=Trlog{y-B) . (2.7) 



operator [13-15 
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Here y is called the boundary cosmological constant, and it parametrizes the moduli space 
of the FZZT brane. Following ||, the simplest such observable is the exact partition 
function of the FZZT brane: 



e W(y) 



(det(y-S)> =Z p , 1 (g)- 1 J dAdB e-^ (v ^ A+1) - AB) det(y - B) . (2A 



Here the integral over A and B is defined in the same way as in ( |2.6| ), as discussed above. 
For the rest of this section, we will focus mainly on the single determinant expectation 
value ( |2.8| ), since this captures the essentials of the matrix model analysis. In section 2.4 
we will also extend our analysis to the correlator of multiple determinants. 

The expectation value of det(y — B) can be computed efficiently using the method 
of orthogonal polynomials (for a review and original references, see e.g. [JT^JT^]). The bi- 
orthogonal polynomials are defined with the normalization P m (a) = a m + . . . and Q n (b) = 
b n + . . ., and satisfy the orthogonality relation 

J dadbe-^ v ^ a+1) - ab) P m (a)Q n (b) = h m 5 m , n , (2.9) 



where again, the integration over a and b is defined as in (|2.6|) . For the orthogonality 



relation (|2.9| ), it is possible to write down explicit, closed-form expressions for P n and Q n : 

Pja) = a n 



Q n (b)= (-g°-y eW'+V-W 



(2.10) 

z=0 



It is straightforward to verify that ( |2.1Up satisfies ( |2.9[ ). 

In terms of these orthogonal polynomials, the expectation value of the determinant is 
simply 



(det(y - B)) = Q N (y) = {~9^\ ^ e^ z+l) ~ yz) 



(2.11) 

z=0 



Now it only remains to extract its double-scaling limit. 

Before proceeding to this, it is instructive to re-derive the result ( |2.11[ ) using auxiliary 
fermions, as was done for the (2, 1) model in ||. Let us introduce fermionic variables Xi 
and xt where i runs from 1 to TV, and write 

det(y-S) = J dxdx^e x ^ y - B)x . (2.12) 



Then, we can integrate out the matrices A and B as follows: 



(det(y-S)} = Z M (g)- 1 J dAdB d X d X ^ e -^ { v A A + i } -AB )+ ^ {y -B )x 

= Z (Pil) (g)- 1 (2ng) N2 J dAd X d X ] e~ L a T ^ A+1)+ ^ H{A + g X X^) (2.13) 



In the last line we have used the fact that Tr(xX^) k = ~{x [ x) k an d the fact that the 
partition function evaluates to 

z P M = (^g) N2 ■ (2.14) 

The next step is to introduce auxiliary parameters s and z and write 

(det(y-B))= /— / d X d X * e^ z +V-y^eT(*+^ 
J 2ng J 



2ng 1 ' 



(2-15) 



By integrating z by parts, the factor of e lsz ^ 9 {is) N can be converted to e lsz ^ 9 {—gd z ) N 
acting on the rest of the integrand. The integral over s will then give rise to a delta function 
constraint S(z), and integrating over z reproduces the answer ( [2.11| ) for the expectation 
value of the determinant. 

The derivation in terms of the fermions leads to the following physical interpretation. 
The matrices A and B describe the open strings between iV condensed ZZ branes. These 
provide the open-string dual to the (p, 1) closed-string background, as shown (for c = 1) in 
|T^,|19|. Introducing the FZZT brane leads to X an d X ) which represent the open strings 
(including orientation) stretched between the FZZT brane and the condensed ZZ branes. 
Finally s and z represent the degrees of freedom for the effective theory on the FZZT brane, 
obtained after integrating out all the background degrees of freedom. It is interesting that 
s and z seem to be conjugate to one another (as do A and B). In the next subsection, 
we will take the double-scaling limit of the s, z integral ( |2.15|) , and we will show that it 
reduces to the generalized Airy function. 



2.2. Double- scaling Limit 



In order to take the double-scaling limit, it is convenient to write (|2.11| ) in yet another 
integral form 

(det(y - B)) = (-g) N N\^- I ^(v p ( z +D-y Z )-Nio gz ^ {2 W) 

2m J z 

where the z contour integral picks up the residue of the pole at z = 0. The integrand 
contains p saddle points in the complex z plane, located at the solutions of 

(z + l) p + yz + gN -1 = . (2.17) 

The saddles collide when y — > and gN = 1, and this is the critical behavior that gives 
rise to the double-scaling limit of the (p, 1) model [l2|. To extract the behavior of ( |2.16| ) 
in the double-scaling limit, we set 

9 = jf, (2-18) 

and scale 

N = e -(p+i) i z = -l + ez, y = e p y, (2.19) 

while sending e to zero. We also deform the small contour around z = to a contour C 
along rays emanating from z = — 1: 

f -e" i7r /2(p+i) t? t < o 

C : z = { (2.20) 
[ e W2(p+i) t5 t > o . 

To see that this contour deformation is valid, notice that the integrand of (2.16) goes to 
zero exponentially fast for — 2 (p+i) < &rgz < 2(p+i) an< ^ 1^1 ~~ ¥ 00 ■ Thus, we can close the 
contour (|2.20| ) at infinity without changing the value of the integral. This closed contour 



encircles z = 0, and since the integrand of (|2.16| ) has no other poles in the finite z plane, the 



contour simply picks up the residue at z = 0. This shows that (|2.20|) is indeed equivalent 
to the original contour around z = 0. 
In the small e limit, ( |2.16| ) becomes 

(det(y -B))= Af(y) J ^ = Af(y)Ai p (y) , (2.21) 

where we have defined 

Ai p (y)= / —e — ~y z , (2.22) 



and M{y) = ^/2nNN~ 1 ^ p+1 ' > e~ N ^ 1 ~ yS> is a non-universal overall factor. Removing this 
factor leads to the final answer for the exact partition function i/j(y) of the FZZT brane in 
the (p, 1) background: 

ij>(y) = lim AT( 2/ )- 1 (det(y - B)) = Ai p {y) . (2.23) 

N — >oo 

This is the generalization of the result of || to higher p. 

In order to avoid cluttering the equations, we will henceforth drop the tildes from the 
double-scaled variables y and z, 

y^y, z^z . (2.24) 

Since we will focus almost entirely on continuum quantities from this point on, we hope 
this notational change will not confuse the reader. 

One recognizes that Ai p (y) satisfies the differential equation 

— J Ai p {y) - yAi p {y) = . (2.25) 

Different choices of contours lead to different linear combinations of the homogeneous 
solutions of (|2.25|) , but the particular contour ( |2.20| ) is special for the following reason. 



It gives rise to the unique solution of ( 2.25Q which is real for real y and decays without 



oscillating at large positive y. (For p = 2, this function is precisely the Airy function Ai(y).) 
It is highly nontrivial that for all p, this physical solution of ( |2.25|) emerges unambiguously 
from the double-scaling limit of the two-matrix model. This is one of the main results of 
our paper. 

By taking y to be large and positive, one recovers the asymptotic behavior 

log^G/) « -p-^y {p+1)/p - ^ lo %y ( 2 - 26 ) 

which we recognize as the large y asymptotics of the FZZT disk and annulus amplitudes 
. Below, we will see that this semiclassical approximation is valid not only for large, 
positive y, but for all large y except y — > — oo. 

Let us also mention that, after taking the double-scaling limit, the contour of inte- 
gration defining ip(y) can be deformed arbitrarily, so long as the contour does not pass 
through "mountains" at infinity. By "mountains," we mean asymptotic regions where the 
integrand grows without bound or decays slower than \z\ . In practice, this means that 
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the deformed contour must asymptote to rays — e zd t and e l9 t as t — > — oo and £ — >■ +00, 
respectively, with 9 satisfying 

n <8< ttt— — — — - . (2.27) 



C: ~t_ v+lt ' (2.28) 



2(p+l) 2(p+l) 

The integral will converge fastest along the path of steepest descent C, which corresponds 

to 6 = i.e. 

p+i ' 

CP+ 1 1, t — ► +00 

Note that the steepest descent contour is at the same time the stationary phase contour. 
In the next subsection, we will use this contour to study the role of Stokes' phenomenon 
in the generalized Airy integral fl2.22f) . 

2.3. Asymptotic expansion and Stokes phenomenon 



Using the explicit integral representation ( 2.22|) and the saddle-point approximation 



it is straightforward to analyze the large \y\ asymptotics of the generalized Airy functions. 
Just as in the case of p = 2, we expect to find a rich structure of Stokes lines and anti-Stokes 
lines. (A succinct account of Stokes' phenomenon can be found in [p5[ .) 
The integral ( |2.22| ) has p saddle points at 

z = z n (y) = \y\hi {e+2irn) , n = 0,...,p-l , (2.29) 

where y = \y\e ld . The value of the argument in the exponential at n-th saddle point is 
simply 

S n = ^z n (yr+ 1 -yz n (y) = -^-\y\ E ^e^ {9+2 ™ ) . (2.30) 
p+1 p + 1 



From the form of the contour (|2.20| ) , it is not difficult to see that for large positive y only 
the n = saddle contributes to the integral. The other saddles are inaccessible due to 
large "mountains" in the integrand of ( |2.22|) at infinity. 

As one varies y in the complex plane, however, other saddles can begin to contribute. 
This happens when the steepest descent contour C collides with the other saddle points. 
The set of points in the complex y plane where this happens form what are called the 
Stokes lines. A necessary condition for this to happen is for two saddle points z n and z m 
to satisfy 

ImS^ = ImS m (2-31) 
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since this means that the two saddles can be connected by a stationary phase contour. 
( |2.31|) is not a sufficient condition for Stokes phenomenon to occur, since we are only in- 
terested the particular stationary phase contour C described in (|2.28|) . In general, though, 



we expect from ( |2.31| ) that the first Stokes line closest to the positive, real y axis occurs at 

(p + 2)tt 



6 



2(p + l) 



(2.32) 



In addition to new saddle points contributing and disappearing, one should allow for 
the possibility of a saddle other than n = dominating. This would happen only when 



Re S n = Re So 



(2.33) 



for some n. Using ( 2.30 ), this becomes 

'p + 1 



cos 



p 



(9 + 2rm) = cos 



p+1 
V 



(2.34) 



whose only solution is 9 = n. The locus of points in the complex y plane where the 
dominant saddle gets exchanged with another saddle are known as the anti-Stokes lines. 
For the generalized Airy functions, the anti-Stokes lines are found to always lie along 
6 = tx. In other words, the asymptotic form of the generalized Airy function is described 
everywhere by the analytic continuation of the n = saddle point, except along the 
negative y axis. Physically, this means that the semiclassical approximation, described in 
( |2.26|) , is valid for all large y except y — > — oo. 

It is useful to illustrate this behavior through a concrete example. Consider the 
contours of fixed phase for p = 5. For small 9, only the dominant n = saddle lies along 
the contour of steepest descent. As one increases 9 from zero to 9 = 77r/12, this contour 
collides with the n = 4 contour. So for 9 > 7n/12, both saddles n = and n = 4 contribute. 
Continuing along in this way, one finds multiple saddles contributing and disappearing, 
summarized in the table below. 



9 = arg(y) 


Contributing Saddle 


< 6 < ^tt 





iiyr < 9 < 7T 


0,4 
0,3,4 
0,2,4 
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A similar result can be obtained for — n < 9 < 0. Finally, at 9 = n, the n = 4 saddle 
replaces the n = saddle as the dominant saddle, while at 9 = — 7r, the n = 1 saddle takes 
over as the dominant saddle. The structure of the saddle points and their interplay with 
the steepest descent contour is illustrated in figure 1. The same analysis can be repeated 



for general p. (The p = 3 case in particular was studied in |2T| in a rather different context.) 



Ott/12 



7 71/12 



9 tt/12 



11 tt/12 



12 tt/12 








Fig. 1: Contours of steepest descent in evaluating the integral expression ( |2.22[ ) for 
the p = 5 generalized Airy function for various \y\ 3> 1 are illustrated in thick (red) 
line. "0,1,2,3,4" label the saddle points. The phases of y are chosen to coincide 
with the Stoke's lines. 



What is interesting about the structure of Stokes and anti-Stokes lines for generalized 
Airy functions is that while the identity of the dominant saddle is different between 9 = n 
and 9 = —n, the actual function is smooth along the anti-Stokes line. The exact answer 
has no branch cut along the negative real axis of the y plane, despite the fact that the 
leading perturbative amplitude ( |2.26| ) is a p-sheeted cover of the y plane. Just as in the 
case of the p = 2 model, the exact FZZT amplitude is an entire function of y, and the 
unphysical sheets of the FZZT moduli space are eliminated by non-perturbative effects. 
However, as we see from the table above, the unphysical sheets still have a physical effect: 
they contribute small, non-perturbative corrections to the exact answer. 

There is one sense in which the p > 2 case is different from the p = 2 case. The 
contribution from the dominant saddle ( |2.30|) arbitrarily close to the negative real axis 

So^-^l^l^ ( 2 - 35 ) 
p+ 1 

has positive real part for p > 2. So Ai p (y) is oscillatory, but it grows exponentially in 
magnitude as y — > — oo. Contrast this with the p = 2 case, where it oscillates with an 
amplitude which decays as a power of y. This difference between p > 2 and p = 2 has the 



following matrix model interpretation. As was discussed in || (see also section 9.3 of ||22|| ), 
the oscillatory power-law decay of Ai(y) as y — » — oo is directly related to the fact that the 
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p = 2 model can be described by the algebraic eigenvalue density of a one-matrix model. 
Applying the same reasoning to p > 2 shows that these models cannot be described by 
an algebraic density of a single eigenvalue - attempting to construct such an eigenvalue 
density from tp(y) = Ai p (y) would lead to a density that grows exponentially as y — > — oo. 
Of course, this is the expected result since, as we saw above, these models are dual to 
two-matrix models which cannot be reduced to one-matrix models. 

2.4. Multiple FZZT 

Finally, let us show how to extend the preceding analysis of the single determinant 
expectation value to the case of multiple determinants. Consider a stack of m FZZT 
branes, whose boundary cosmological constants are described by an m x m matrix Y 
with eigenvalues y%, . . . , y m . The multiple FZZT partition function in the matrix model is 
formulated as a double-scaling limit of 

(det(Y ®I N -I m ® B)) = ^iQ^lM , (2 .36) 

where i and j take on values 1, . . . , m and A(y) is the Vandermonde determinant. The 
derivation of ( |2.36| ) can be found in e.g. [I7|] . 

To evaluate this expression in the large N limit, let us note that according to (j2.18|) , 
a shift in N is equivalent to a shift in g. Now, g is the coupling of the lowest dimension 
operator, and for the (p, 1) models, a shift in this coupling is equivalent to a shift in y. 
Therefore, (|2.36|) becomes in the double-scaling limit 



/ \ det ' ' ib ( u ' 1 

(det(Y ® I N - I m ® 5)) -> if, m (Y) = tJ ^ y) ■ (2-37) 

Here the finite and continuum Y are related exactly as in ( |2.19|) , and we have dropped 
irrelevant factors just as in ( |2.23| ). This formula for the multi-FZZT correlator is a special 
case of the general formula (applicable for all (p, q) models) in 0. 

Using (|2.23|) along with the integral representation (|2.22|) , we can turn (ET37| ) into 



l> m (Y) = AipiY) , (2.38) 
where Ai p (Y) is the generalized matrix Airy function 

Ai p (Y) = A(y)- 1 J d Zl A(z)ei^- mz * = J dZe Tr (^~ YZ ) . (2.39) 

To establish the last equality, we used an identity from [p3| , |2^| . Eq. (|2.38|) shows explicitly 
how the effective theory on m FZZT branes in the (p, 1) background is described by the 
m x m generalized Kontsevich integral. 
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3. Open/closed string correspondence 

The main result of the previous section can be summarized as the statement regarding 
the expectation value of the exponentiated macroscopic loop operator 

(e w ^)=^{y)=Ai p {y) (3.1) 

in the double-scaling limit. These macroscopic loop operators, when parameterized by the 
length I of the boundary on the world sheet, can be decomposed into local closed- string 
operators in terms of their scaling [ ]TB| - |To| , 

pk/p 

fc>i r \~^) 

where Clp IS ct dimensionless normalization whose precise value will be fixed in appendix B. 
The I and the y parametrization of the macroscopic loop operators are related by Laplace 
transform 

W(y) = I je- e yW(£) ~ £ V^a^a, . (3.3) 

J k>l 

A similar relation between boundaries of fixed length and local operators (as well as a 
Laplace transformed version) also appears in section 3.1 of ||10|1 . One should keep in 
mind, however, that the expansion of macroscopic loop amplitudes in terms of microscopic 
operator correlation functions is subtle, in that there are divergent contributions as I = or 
y — > oo in the disk and the annulus amplitudes [[14] . These can be removed by introducing 
the normalization factor 

^, x 1 P_ y (P+i)/P_ (P- 1 ) logif fn 

C(y) = —j=e p+ lV 2 p gy , (3.4) 

which precisely cancels the divergences from the disk and the annulus amplitudes. One 
can then interpret the result of the previous section as the statement 

k/p a p , k a k ) ) = C(y)-V(y) = C{y)- x AUy) , (3.5) 




which relates the normalized FZZT brane partition function to the generating function of 
closed-string correlators in (p, 1) topological string theory. 



The relation (|3.5| ) is the simplest statement of open/closed string duality - it says that 



the partition function of the FZZT brane is equivalent to a certain shift in the closed-string 
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background. This relation becomes much richer if one considers multiple FZZT branes, as 
in section 2.4. The left hand side of (|3.5|) simply generalizes to 



exp 



while the right hand side of ( |3.5| ) becomes 

Z m (Y) = C(Y)- 1 Ai p (Y) . (3.7) 
Here Ai p (Y) is the matrix Airy function defined in ( |2.39| ), and 

C(Y)= lim Ai p (Y) = dZ e { 1+p + a2^=o y zy z ) (3.8) 

with the contour of Z integration defined as in (|2.39|) . Combining all these ingredients, we 
arrive at the statement of open/closed string duality 



exp 



J2tk°k ) = S m (F), t k = lTrY- k / p a p , k , (3.9) 



k 



which relates the generating function of topological closed string amplitudes to the gener- 
alized Kontsevich integral. 

Let us note that another way of writing ( 3 -71) , which may be more familiar to some, 
comes from substituting Y = X p 

ZP+1--XP+ 1 



Tr ^ ' ' -(Z-X)X p 

Sm (X) = JdZe ' — . (3.10) 

J dZe ^i=0 2 

This form of the generalized Kontsevich integral is common in the literature (see e.g. ||16| ), 
modulo trivial rescalings and shifts of Z and X. In the worldsheet (Liouville) description of 
these models, Y corresponds to the boundary cosmological constants of the FZZT branes, 
and X corresponds to the dual boundary cosmological constants. 



4. Discussion 

In this article, we computed the exact partition function of (multiple) FZZT branes in 
the (p, 1) topological background, using the two-matrix model in the double- scaling limit. 
We found that these partition functions are given by the generalized Kontsevich integral. 
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By relating this to a specific insertion of closed string operators, we were able to formulate 
a precise open/closed string duality for (p, 1) topological gravity with matter. 

In principle, the double-scaled matrix model can provide a framework for understand- 
ing open/closed string duality in the most general minimal string theory. However, one 
should keep in mind that our analysis was aided by various simplifications that occur in the 
(p, 1) models. The most important simplification came in the calculation of the FZZT par- 
tition function at finite N. In general, the FZZT partition function is the scaling limit of an 
orthogonal polynomial of the dual matrix model. For the (p, 1) models, these orthogonal 
polynomials were sufficiently simple as to admit an elementary, closed-form representation. 
Furthermore, this representation allowed the scaling limit to be taken explicitly, which led 
directly to the generalized Airy function. 

It would be interesting to push the matrix model analysis of open/closed string duality 
to more general (p, q) models which are non-perturbatively well defined. Of course, the 
orthogonal polynomials for the general (p, q) models are more complicated. Nonetheless, 
they are known to satisfy a recursion relation and can be generated in a finite number 
of steps. So there is still hope that one may be able to formulate concretely open/closed 
string duality for general (p, q) . 

Another open problem of immediate interest is to rederive the results in this paper us- 
ing open-string field theory, as was done for (2, 1) in p. (See also [p5[ for recent discussion 
on this issue.) This would presumably shed more light both on open/closed string duality 
in these models, and on the inner workings of open-string field theory itself. It would also 
confirm the identification of the generalized Kontsevich integral with the effective theory 
of open strings between FZZT branes. 

The primary motivation for studying explicit realizations of open/closed string duality 
in toy systems like (p, 1) is to provide new insights that can be extended to richer dynamical 
systems. Let us briefly mention one possible lesson that could be learned from our work. 
So far, all attempts to formulate purely closed-string observables from open-string field 
theory in general [26 , 2j| (see also the interesting recent work of [28]), following the work 
of fl2T||n|| , have had difficulty in removing the boundary of the world sheet. By contrast, 
here we encountered no such difficulties in formulating the open/closed string duality of 
the (p, 1) model: purely closed-string amplitudes were easily obtained from the generalized 
Kontsevich matrix integral, and the latter presumably represents a reduction of open-string 
field theory along the lines of ||. Nonetheless, there was one seemingly unnatural step 
where we removed the divergent contribution of the disk and the annulus amplitudes by 
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hand. Similar truncations of a few terms on the open string side also arose in the context 
of Gopakumar-Vafa correspondence ||32|| . It would be interesting if such truncations are 
significant and have ramifications for open/closed string dualities in general. 
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Appendix A. Generalizing to other potentials 

In this article we studied the two matrix model with a single multi-critical potential. 
However, it is very easy to generalize this discussion to superpositions of multi-critical 
potentials. This will allow us to use the two-matrix model to study the integrable flows 
between multi-critical points with different values of p, an aspect of minimal string theory 
that could not be studied in |J where the one-matrix model (which always describes p = 2) 
was used. 

So let us consider the matrix integral 




(A.l) 



where the potential is determined using the formula ( |2.2|) of |T^] 



z 




(A.2) 



k=0 



Then the expectation value of the determinant becomes 



(det(y - B)) 



2ni J z 



V k (z+l))-yz) 



-N log 2 



(A.3) 
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One can then take the double-scaling limit 

p-2 

N = e~ p -\ y = e p y, z = -l + ez, gN=l + Y,^ P ~ k s k (A.4) 

k=0 



Then as e — > 0, (|A.3| ) becomes 



e- Ny (det(y-B)) ~ — / dze^+ELN^-S* = ^(y) . (A.5) 

where we have dropped numerical factors which do not depend on y, and the contour are 
defined exactly as before. 

One can show explicitly that the FZZT partition function ( |A.5| ) is the Baker- Akhiezer 
function of the KP hierarchy with Lax pair 



Q = d p + J2 s kd k , P = Q 1 l P = d, (A.6) 



>i/p 

k=0 

where d = d So = —d y . These clearly satisfy 



[P, Q] = 1 (A.7) 
and one can easily check that if) satisfies the equations defining the Baker- Akhiezer function: 

QtP = ffl>, PtP = -d y tp . (A.8) 

Note that we can convert this to a differential equation for ip(y) as a function of y, 

p-i 

i> ip \y) + W fc) (!/) = (-i) p yV(y) (A.9) 

which is a further generalization of the Airy equation. 

Appendix B. Generalized Kontsevich integral and intersection numbers 

Here we will fix the normalization constants a Py k that arose in the decomposition ( |3.2| ) 
of the macroscopic loop operator into closed-string scaling operators. Our method will be 
to match H m (X) given in ( |3.10| ) with the generating function 



n . /n 



r(t) = exp(F) = exp ( £ < J] <-> J] ^ ] (B.l) 

{n,m} {n,m} {n,m} ' 
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of intersection numbers associated to (p, 1) topological string theory. These were defined 
in |33[] as the integral 

S 

= ^[- tlci(£i) ni cn(V) (B.2) 

P J M g , s (p,m) i=1 

over the stable compactification of a covering jM SjS (p, m) of the moduli space A4 g , s of a 
Riemann surface E of genus g with s distinct marked points Xi, i = 1, . . . , s. Here a"o,m, 
to = 0, . . . , p — 2 are the primary matter fields and n = 0, . . . labels the n-th gravitational 
descendant a n>m . V is a vector bundle over .M 9)S (p, to) with fiber H°(C,1Z) of dimension 
D = (g - 1)(1 - |) + ^ ^ and c D (V) is its top Chern class. The line bundle "ft is 
defined up to isomorphism by 1Z®p = K p ~ 1 &) s i=1 0(xi)® mi , where K is the canonical 
bundle of £ and 0{x{) are line bundles whose sections have at most simple poles at Xi. 
The covering M. g ^ s {j>, m) depends on data m,p as indicated by the notation. ci(£j) are 
the usual Mumford-Morita-Miller classes. On dimensional grounds the integral above is 
zero unless 

X> + y - l) = (s - (l - £) ) (9- 1) • (B.3) 



A mathematically precise construction of cd(V) was formulated in p3,B3]. The authors 



|34[1 furhter confirmed the explicit calculation of intersection numbers in |33]]. A suitable 



compactification of A4 9iS (p, to) has also been discussed in |36],|37]. 



As was conjectured by JT^,|TTj] and proven for the case of p = 2 in [|10|j , ( p.lj ) is the 
tau function of the p reduced KP hierarchy fulfilling in addition the string equation. The 
geometric intersection numbers ( P-2|) come with natural normalization, and the first few 
were worked out in ]33|]. We shall use ( [B.l|) with this normalization of the intersection 
numbers as the definition of the closed string couplings t n m . 



Following [PBf , one can show that F = lim m ^ 00 log(S m ) fulfills the string equation 
and the p reduced KP hierarchy in the Miura variables tk ~ Tt(X~ k ). These properties 
determine F up to an additive constant. One can choose the normalization of tk so that 



the string equation 



dF 1 ^ OF 



5 Y £ 0,i*0,j + Y Y Wl.m^T (B.4) 



x+j=p—2 n=Om=0 
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and the equations of the p reduced KP hierarchy 

= -c„, fc res (Q n+h ^) , Q = DP - V^^^})^, D = -^pj- 

(B.5) 

take on precisely the form used in [p3| , i.e. with c n> k = (—l) n p n+1 / Ylj=o(iP Jr ^ + -*-)■ This 
can be achieved by setting 

t np+fc+1 = t n>fc = (-;p) fe ~V + ( i P ) +2) J](^ + fc + l)Tr(X-^+ fc+1 )) . (B.6) 



This fixes the normalization constant a P; k in (|3.9| ). 

As a check of our analysis, we have calculated a few of these intersection numbers 
explicitly from the asymptotic expansion of S m (X), and compared them with results in 
the literature. To extract the intersection numbers from S m (X), we follow ||38|| . the only 
difference being that the differential operator D with the defining properties 



J dz ze P+ 1 



p+1 



or (D p — = 0, is given in our parametrization by 



n p-11 Id 

D = + x — — . (B.J 

2p x p px p 1 dx 



The expression for S m (x) in terms of generalized characters of the symmetric group in ||38|| 
gives the exact result for all coefficients involving Tr(x _fc ) up to k(p+ 1) < m. Expanding 
in the corresponding order and using ( |B.6| ) we obtain 



1,2 1 

'/<-•"> - \ 2*0,0*0,1 + T2 fcl '° 



^0=3 — ( -tn nto.l + —ti 

(B.9) 



■ 1 4 1 2 1 2 1 3 1 

+ ( 72*0,1 + 2*0,0*0,1*1,0 + 24*1,0 + g*o, 0*1,1 + ^0,0*2,0 J + 



and 



/1 2 1 2 1 \ / 1 2 2 1 2 1 2 

p=4 ~~ I 2 to '° toA + 2 to, ° t0,2 + 8* 1 ' I + \ 16* ' 1 * ' 2 + 2* ' * - 1 * 1,0 + 2 to '° to ' 2tl ' 

1 2 1 2 1 3 1 1 \ 



(B.10) 



1 In the case of p = 3, the same change of variables was described in [39], 
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for p = 3 and p = 4 respectively.0 In particular, we find agreement with the explicit 



calculations of [33] for the chiral ring correlators o k = ctq k 



(^^2^3) — S kl +k 2 +k 3 ,p-2 

(a kl a k2 a k3 cJk 4 ) = -mm(ki,p - 1 - h) , 
p 



(B.ll) 



from which all chiral ring correlators on the sphere follow by associativity. We also agree 
with (r^o) = (p — l)/24. More generally we find that for p < 10 and k < 3, all predicted 
intersection numbers are positive rational numbers. Further predictions for the intersection 
numbers with homogeneous degree d where a n> k has degree pn + k + 1 are tabulated below 
for p = 3, 4, 5. 



p = 3 


d = 12 


(c r o,i°"i.o)o 




2 
3 


(c r o,o cr o,i <J i,o 





2 


(cr 00 (T2,i)o 




1 


(O"0,00"o 1(71,1 





1 

3 


(C0,0 (7 0,1 <T 2,0 





1 


(o"o,oO"i,oo-i,i 





2 






1 






6 


(o"0,lO-i,i)l 




1 

36 


(o"0,0O" l,OC"2,o) 


1 


1 

6 


(<jQ,l(T2,x)l 




1 

36 


(O"0, 0°"3,o)l 




1 

12 



p = 4 


d = 15 


((JO i O°"2,2)o 


1 


(Co,0 "0,l°"2,l)o 


1 


(Co,0°"0,20"2 > o)o 


1 


(°0,0 0-0,l (J 2,o)o 


1 


(CO.O "!, 0^1,2)0 


2 


(cO^CTo,! CTq, 2^1,2)0 


1 

4 



p = 4 cont. 


d= 15 




p = 5 


d = 6, 12 


( cr o,i fJ i,2)o 


1 
4 




(an 00-0,3)0 


1 


( cr o,o <J i,i)o 


2 




((To,00-0,lO-0,2)o 


1 


( cr o,o cr o,i fJ i,oo"i,i)o 


2 






1 


(o"o,00"0,lO"o ) 20"l i l)o 


1 

/l 

4t 




(Oq, 00-1,3)0 


1 


(Co,l (J 0,20"l,l)o 


1 
2 




(Oo,oO"0,lO-l,2)o 


1 


/ 2 2 \ 
\C0, 0^0,20-1^/0 


2 




(oo,ooo, 20-1,1)0 


1 


/ 2 2 \ 
\CTO ) 00'o > lO' lj0/ )o 


2 




(00,000,10-1,1)0 


1 


( fJ 0,l fJ 0,20"l,o)o 


1 
2 




( a o,o 0-0,30-1,0)0 


1 


(O"o,2)o 


1 

8 




(o"o,oO"o,lOo,20-l,o)o 


1 


(00,00-1,2)1 


1 

48 




(00,10-1,0)0 


1 


(^cU^o)! 


1 
8 




( a o,i ^0,3)0 


1 
5 


(0-0,000, 202,2)1 


1 

96 




(00,10-0,20-0,3)0 


1 

5 


(0-0,10-2,2, )l 


1 

32 




(0-0,2)0 


2 
5 


(o"o, 10-0,20-2,1)1 


1 

24 




(o-0,00"2,o)l 


1 

6 


(0-0,00-1,00-2,0, )l 


1 
4 




(0-0,2 0"1,3, )l 


1 

60 


(O"o,20"2,o)l 


1 

48 




(0-0,30-1,2)1 


1 
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(o-l,o)l 
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48 
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6 
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1 
4 








(0-0,20-1,1)1 


1 

24 








(0-3,2)2 


3 

2560 









For p = 2 our expansion agrees with the one in [38]. 
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Finally, let us emphasize that, despite all the progress in making Witten's conjecture 
mathematically precise [ j34H37|1 , a proof of (P-l| ) is still lacking for p > 2. That is, it remains 
to extend beyond p = 2 Kontsevich's statement that the intersection numbers for (p, 1) 
topological gravity are generated by the tau function of the p reduced KP hierarchy. Such 
a generalization will presumably involve open string field theory techniques in order to 
generate the cell decomposition of the moduli space of spin curves. See [P|,|23[] for a recent 
discussion on this issue. 
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